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ABSTRACT 

A recent observation with the Hipparcos satelhte and some numerical simulations 
imply that the interaction between an oblate halo and a disc is inappropriate for the 
persistence of galactic warps. Then, we have compared the time evolution of galactic 
warps in a prolate halo with that in an oblate halo. The haloes were approximated as 
fixed potentials, while the discs were represented by A^-body particles. We have found 
that the warping in the oblate halo continues to wind up, and finally disappears. 
On the other hand, for the prolate halo model, the precession rate of the outer disc 
increases when the precession of the outer disc recedes from that of the inner disc, 
and vice versa. Consequently, the warping in the prolate halo persisted to the end 
of the simulation by retaining the alignment of the line of nodes of the warped disc. 
Therefore, our results suggest that prolate haloes could sustain galactic warps. The 
physical mechanism of the persistence of warp is discussed on the basis of the torque 
between a halo and a disc and that between the inner and outer regions of the disc. 

Key words: galaxies: haloes - galaxies: kinematics and dynamics - galaxies: structure 
- methods: numerical 
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1 INTRODUCTION 



Many spiral galaxies, including our Galaxy, have warped 
discs which resemble characteristic 'cosmic integral signs'. 
That is, the outer disc lies above the inner disc plane on one 
side, and falls bellow that plane on the other side. Although 
the warping is often seen in neutral hydrogen layers (San- 
cisi 1976; Bosma 1981), it is also observed in stellar discs 
(van der Kruit & Searle 1981; Innanen et al. 1982; Sasaki 
1987). In the Milky Way, the stellar warp has been detected 
not only for young stars (Miyamoto, Yoshizawa & Suzuki 
1988) but for old stars (Porcel, Battaner & Jimenez- Vicente 
1997). In addition, the frequency of warped discs in spiral 
galaxies is sufficiently large that at least half of spirals are 
warped both in H I discs (Bosma 1991) and in optical discs 
(Sanchez-Saavedra, Battaner & Florido 1990; Reshetnikov & 
Combes 1998). These observations imply that warps must 
persist for a long time unless they are repeatedly excited. 
It is true that some galaxies with warped discs (e.g., M31, 
see Innanen et al. 1982) have nearby companions. However, 
there do exist warped galaxies (e.g., NGC 4565, see San- 
cisi 1976) that have no nearby companions being supposed 
to be responsible for the warp in the recent past. In fact, 
Reshetnikov & Combes (1998) have revealed that about 21 



* E-mail: ideta@kusastro.kyoto-u.ac.jp (MI); 

hozunii@sue.shiga-u.ac.jp (SH); tsuchiya@kusastro.kyoto-u.ac.jp 



(TT); takizawa@kusastro.kyoto-u.ac.jp (MT) 



out of 133 isolated galaxies are warped like an integral sign. 
This indicates that warps are not necessarily caused by tidal 
interactions with other galaxies. 

One explanation for isolated warped galaxies is the 
gravitational torque of a halo acting on a disc that is 'mis- 
aligned' to the equatorial plane of the halo. Such a tilted disc 
embedded in a halo has intrinsic spin, so that it precesses 
like a top. Since the precessing rate is a function of radius, 
kinematical warps will wind up and disperse in a short pe- 
riod of time. Once the self-gravity of the disc is taken into 
account, realistic warped configurations emerge in which a 
disc precesses coherently like a solid body inside axisym- 
metric haloes (Sparke 1984; Sparke & Casertano 1988, here- 
after SC; Kuijken 1991). Even if a tilted disc is formed in 
a halo as a different shape from a warped mode, it will be 
finally turned into the mode within a Hubble time (Hofner 
& Sparke 1994). However, according to some numerical sim- 
ulations (Dubinski & Kuijken 1995; Binney, Jiang & Dutta 
1998), the warping in oblate haloes is not retained persis- 
tently, and so, disappears within a few dynamical times. 
Thus, the interaction between an oblate halo and a disc will 
be inappropriate for long-lasting warps. 

Recently, Smart et al. (1998) have extracted warp- 
induced motions in the Milky Way by analysing the data 
obtained with the Hipparcos satellite. Then, they have found 
that the Galactic warp rotates in the same direction as the 
Galaxy. As shown by Nelson & Tremaine (1995), oblate 
haloes lead to the opposite sense of the warp precession to 
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the Galactic rotation, whereas prolate haloes make them ro- 
tate in accordance with each other. In addition, they have 
demonstrated that in some cases, prolate haloes can excite 
warps. Thus, prolate haloes are favourable to the explana- 
tion of Smart et al.'s (1998) finding, if the motions that 
they found are attributed to the interaction of the Galac- 
tic disc with an often assumed massive halo. Cosmologi- 
cal simulations, based on a cold dark matter scenario, have 
also revealed that dark matter haloes surrounding individual 
galaxies are highly triaxial and that the fraction of prolate 
haloes is roughly equal to that of oblate haloes (Dubinski & 
Carlberg 1991). In spite of these circumstances mentioned 
above, prolate haloes have somehow often been ignored in 
previous studies on warps. Therefore, we need to pay atten- 
tion to the warp arising from prolate haloes. 

In this paper, we examine how a warp is developed and 
evolves in prolate haloes in comparison with that in oblate 
haloes. As a first step, we treat the haloes as external fixed 
potentials. In Section 2, we describe the models and the 
numerical method. Results are presented in Section 3. In 
Section 4, we analyse our results and explain them on the 
basis of a simplified model. Conclusions are given in Section 
5. 



2 MODELS AND METHOD 

We study the evolution of self-gravitating discs embedded 
in axisymmetric haloes. As shown by Nelson & Tremaine 
(1995) and by Dubinski & Kuijken (1995), dynamical fric- 
tion between a disc and a halo plays an important role to 
precessing bending modes at least for the inner region of the 
composite system. However, the accurate evaluation of dy- 
namical friction would require a prohibitively huge number 
of particles to represent the halo as well as the disc. Oth- 
erwise, the disc will thicken owning to two-body relaxation 
originating from Poisson fluctuations. In fact, Dubinski & 
Kuijken (1995) reported the vertical disc thickening for a 
self-consistent model with a particle disc, bulge, and halo. 
As a result, warped structures developed in the disc could 
not be distinguished from the background particle distribu- 
tion, which would lead us to an incorrect conclusion about 
the longevity of the warp. Then, we begin with rigid halo 
models, as a first step, to unravel the effects of the halo 
shape on the warp. 

The density distribution of the halo is an axisymmetric 
modification of Hernquist's models (Hernquist 1990) being 
suitable for spherical galaxies and bulges. Then, the halo 
density profile is represented, in cylindrical coordinates, by 

Ph{R,z) = - — 5 -3, (1) 

where Mh is the halo mass, a and c are the radial and vertical 



core radii, respectively, and 
R^ 



2 d'z^ 



(2) 



The cumulative mass profile and potential of the halo are 
written, respectively, by (Binney & Tremaine 1987), 



Mh [R, 



and 



Mh 



(1 + ^ 



(3) 



Table 1. Parameters for Af-body simulations. 
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Figure 1. Circular speed showing the contribution from the disc 
and halo to the total for the oblate (top) and prolate (bottom) 
halo models. 



$h {R, z) = 



GMh 



du 
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where 



R^ 



+ 



a^ -I- M c^ -f It 



(4) 



(5) 



As a realistic disc model, though a bulge component 
is not included, we adopt an exponential density profile in 
the radial direction (Freeman 1970) and an isothermal sheet 
approximation in the vertical direction (Spitzer 1942) given 
by 



Pd{R,z) 



Md 



■ exp 



V Ra)^^^ Ud/ 



(6) 



where Md is the disc mass, R^ is the disc scale- length, and 
Zd is the disc scale-height. The discs are truncated radi- 
ally at 15 -Rd and vertically at 2 z^- Following Hernquist's 
(1993) approach, we approximate the velocity distribution of 
disc particles using moments of the coUisionless Boltzmann 
equation; the velocities are sampled from Gaussian distri- 
butions with means and dispersions derived from the Jeans 
equations. The Toomre Q parameter (Toomre 1964) used to 
normalize the radial velocity dispersion is set to be 1.5 at 
the solar radius, Rq = (8.5/3.5) R^- The Q profile varies 
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with radius and has a minimum at approximately R — 2 R^ , 
where the minimum Q value is about 1.47. To avoid compli- 
cations due to an extra component such as a bar, this rather 
large Q distribution is chosen so that the bar instability will 
not occur in the disc. 

We have ensured that the disc-halo system thus con- 
structed is really in equilibrium if the disc is initially placed 
in the equatorial plane of the halo, because the density pro- 
file of the disc did not change significantly over several or- 
bital times. However, the disc does not necessarily form in 
the equatorial plane of the halo. In fact, Katz & Gunn (1991) 
showed that the disc was misaligned to the symmetry plane 
of the halo at an angle of typically 30 degrees at its birth. 
Therefore, in our simulations the disc is initially tilted by 30 
degrees with respect to the equatorial plane of the halo. 

We employ a system of units such that the gravitational 
constant G = 1, the disc mass Md ~ 1, and the exponential 
scale- length _Rd = 1. The orbital period at the half-mass 
radius of the exponential disc, R ~ 1.7 -Rd, is 13.4 in our 
system of units. If these units are scaled to physical values 
appropriate for the Milky Way, i.e., Rd — 3.5 kpc and Md = 
5.6 X 10^" M0, unit time and velocity are 1.31 x 10^ yr and 
262 km s~^, respectively. 

The disc is represented by 100 000 particles of equal 
mass. We show the parameters of our models in Table hi 
and the rotation curves of each model in Fig. hi. The halo 
mass is determined so that the disc and halo masses within 
3 Rd are equal to each other. 

The simulations are run with a hierarchical tree algo- 
rithm (Barnes & Hut 1986) using the GRAPE-4, a special- 
purpose computer for gravitationally interacting particles 
(Sugimoto et al. 1990; Makino et al. 1997). We adopt an 
opening angle criterion, 6 — 0.75. Only monopole terms are 
included in the tree code. The equations of motion are inte- 
grated with a fixed time-step. At = 0.1, using a time-centred 
leapfrog method. The Plummer softening length is 0.04 _Rd, 
or in other words, 0.2 Zd- 



3 RESULTS 

We stopped the simulations at t = 400. This time corre- 
sponds to about 30 orbital periods at the half-mass radius 
of the disc. No bar instability was found in the discs. In 
either simulation, the total energy was conserved to better 
than 0.2 per cent. 

We measured the inclination and the longitude of as- 
cending nodes of the disc relative to the equatorial plane of 
the halo by calculating the principal moments of inertia for 
the bound particles. In Fig. H, we then show the evolving 
density profiles from an edge-on view of the discs in the pre- 
cessing frame. In this frame, an observer is always on that 
line of nodes of the disc which is calculated for the parti- 
cles within the half-mass radius of the disc, R ~ 1.7 Rd- 
The precession periods of the discs in the oblate and prolate 
haloes evaluated with least-squares fits are Toh = 266 and 
Tpr = 306, respectively. These values are in excellent agree- 
ment with those predicted by linear theory [see equation 
(21) of SC] which gives Tob = 266 and Tp, = 308. 

Hofner & Sparke (1995) showed that warped configura- 
tions are developed in oblate haloes. We further find from 
Fig. y that such configurations appear in the prolate halo as 





^ ^^ 
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Figure 2. Time evolution of the projected density distributions 
of the disc in the precessing frame. The top (bottom) two rows 
correspond to the oblate (prolate) halo model. The time is shown 
at the right-up corner of each frame. 

well as in the oblate one. In addition. Fig. g demonstrates 
that the shape of warp depends on that of halo: using the 
terminology in SC, the type I warp that bends upward away 
from the symmetric plane of the halo is developed in the 
oblate halo, while the type II warp that bends down toward 
the symmetric plane of the halo is developed in the prolate 
one. The warp in the oblate halo decayed and disappeared 
almost completely by the end of the simulation, while the 
warp in the prolate one persisted to the end. 

To see the differential precession of the disc, we divided 
the distance between the centre and 5 -Rd evenly into 10 
annuli, and calculated the inclination and azimuthal angles 
of each annulus which contains at least 2 000 particles. Fig. H 
shows, on the polar diagram, the line of ascending nodes of 
each annulus with respect to the equatorial plane of the halo. 
There are two differences between the oblate and prolate 
halo models. One is the sense of the precession; for the oblate 
halo the warp rotates in the opposite direction to the disc 
rotation, while for the prolate one it rotates in the same 
direction. The other is the behaviour of the winding of the 
warp; for the oblate halo the warp winds up tightly with 
time, while for the prolate one the longitude of each annulus 
is kept almost aligned. 

We find from Fig. H that the inner region of the disc, 
-R ;$ 3 7?d, precesses at almost a constant rate independent 
of radius both in the prolate and oblate haloes. Since the disc 
mass is equal to the halo mass within 3 -Rd for our models, 
the self-gravity of the disc is dominant as compared to that 
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OBLATE (T=240 - T = 400) 



PROLATE (T=240 - T = 400) 



Figure 3. Time evolution of the line of ascending nodes of discs 
on polar diagram. The left and right columns correspond to the 
oblate and prolate halo models, respectively. Concentric dotted 
lines are drawn at a radial interval of one sale length, with a solid 
line at five scale lengths. The disc rotation is counterclockwise. 



of the halo within such radius. We thus understand that this 
behaviour arises from the predominance of the self-gravity of 
the disc over that of the halo, as shown by Lovelace (1998). 
On the other hand, the outer disc (i? ^ 3 Ra) precesses at 
a different rate from radius to radius. 

To make clear the difference in precession at large radii 
between the oblate and prolate halo models, we present in 
Fig. W the time evolution of the longitude (top row) and 
that of the inclination angle (bottom row) with respect to 
the equatorial plane of the halo for the outer and inner discs 
which correspond to the annulus between 4.5 Ri and 5.0 Rd 
and to that between 1.5 7?d and 2.0 7?d, respectively. Fig 
shows that the precession rate of the outer disc increases 
with decreasing inclination angle, and vice versa. For the 
prolate halo, at the beginning of the simulation, the incli- 
nation angle decreased, and the precession rate increased. 
In the subsequent evolution, the longitude of the outer disc 
passed through that of the inner disc as the inclination de- 
creased. After the passage of the longitude of the outer disc, 
the inclination increased, and the precession rate decreased. 
For the oblate halo, on the other hand, the inclination in- 
creased initially, and the precession rate decreased. The dif- 
ference in longitude between the inner and outer discs be- 
came larger with increasing inclination, and so, the warp 
wound up with time as seen in Fig. g. 

The difference in inclination angle between the inner 
and outer discs in the oblate halo is larger than that in the 
prolate one at the end of the simulation. However, for the 
oblate halo, the warp disappeared as seen in Fig. bl because 
the azimuthal angle of the inner disc was different from that 
of the outer disc. 
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Figure 4. Time evolution of the longitude of ascending node 
(top row) and the inclination (bottom row) of discs. The solid 
and dashed lines show the outer and inner discs, respectively; the 
inner disc corresponds to a radial annulus between 1.5 R^ and 
2 iJd, and the outer disc to that between 4.5 R^ and 5 R^. 



4 PHYSICAL INTERPRETATION 

As found in the previous section, the prolate halo is plausible 
for the maintenance of galactic warps in that the winding 
problem is avoided. However, there remains a question as to 
what makes the difference in evolution of warps between the 
oblate and prolate halo models. 

Since Lovelace (1998) has shown that the self-gravity 
of the disc can synchronize the precession rate in the in- 
ner region, we need to explain the different behaviour of the 
warp in the outer region between the oblate and prolate halo 
models. Then, we simplify the A''-body models used in our 
simulations and construct a three-component system con- 
sisting of an outer disc, an inner disc, and an axisymmetric 
halo in order to pay special attention to the torque between 
the inner and outer discs in addition to that from the halo. 
Here, we approximate an outer disc as a ring. 



4.1 Simple model 

In this subsection, we solve the equation of motion for the 
outer ring in order to examine whether the results found 
in the A^'-body simulations can be reproduced. For this pur- 
pose, we consider the axisymmetric Binney (1981) potentials 
as models of the halo and the inner disc, 



$h = iK\ln(i?c.h + i?' + ^ ), 

$d = \vI^\ArU + r^^^-^\, 



(7) 
(8) 



where 14 is the asymptotic circular velocity, R^ is the core 
radius, and q is the potential flattening with the subscripts 
'h' and 'd' denoting the halo and the disc, respectively. The 
potentials of the halo and the inner disc are fixed. The inner 
disc is tilted by an angle of 30 degrees relative to the sym- 
metric plane of the halo and given a constant precession rate 
that is calculated from linear theory (SC). The dynamics of 
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Table 2. Parameters for three-component models. 
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Figure 5. Same as Fig. H, but for three-component models. 



the outer ring are solved on the basis of Euler's equation of 
motion for a rigid body (Goldstein 1980). 

The radius and angular speed of the outer ring, Q, are 
set to be 5.0 and 0.15, respectively. The parameters of this 
model are summarized in Table H. The values of the poten- 
tial flattening, Qh and Qd, are adjusted to those which are 
evaluated from the halo and disc models at 5 Rd employed 
in the A'^-body simulations (see Table hi). We determine the 
values of 14, d and 14, d in the same manner. 

Fig. ra presents the time evolution of the longitude and 
inclination of the outer ring, which corresponds to Fig. H 
in the A''-body simulations. We can see that the precession 
rate of the outer ring increases with decreasing inclination 
angle with respect to the equatorial plane of the halo, and 
vice versa. Moreover, for the prolate halo the inclination 
angle decreases at the beginning, while for the oblate one it 
increases. Thus, for the prolate halo, the precession of the 
outer ring can pass through that of the inner disc. After the 
passage, the inclination angle increases and the precession 
rate decreases. This means that the longitude of the inner 
disc and the outer ring remains almost aligned. Therefore, 
the behaviour of the outer ring is the same as that seen 
in the A'^-body simulations. Since we have found that the 
simple model can reproduce the main properties of the A''- 
body simulations, we can rely on this model to figure out 
the physical mechanism of warps in the oblate and prolate 
haloes as described below. 



4.2 The precession 

The precession of the outer ring is caused by the torque due 
to the halo and the inner disc. We take a coordinate system 
in which the z-axis is along the symmetry axis of the halo, 
and the a;- axis is along the line of nodes of the outer ring 




Figure 6. Geometry for the outer disc. The line of node is along 
the z-axis, and the symmetry axis of the halo is along the z-axis. 



(where the outer ring intersects with the z = plane). The 
xy-plane is in the inertial frame. The geometry is shown in 

Fig- I- 

The x-component of the torque exerted by the halo on 
the outer ring is given by 

i'2tt 

Th,^ = - / rd(PX {Fi,,yZ ~ Fh,,y) As-—, (9) 

Jo 27rr- 

where A is the line density of the outer ring, and m is the 
mass of the outer ring. The y- and z-components, i4,y and 
Fh,z, of the force due to the halo per unit mass are, respec- 
tively, written by 



"h,y 



= -^hV, -Fh.z = 



2 



(10) 



where Q,h and /ih are the halo contributions of the orbital 
and vertical frequencies, respectively. The position vector of 
a point on the outer ring is 



R 



(r cos 0, r cos i sm 0, r sm i sm i 



< < 27r, 



(11) 



where i is the inclination angle of the outer ring with respect 
to the equatorial plane of the halo. Then, on the assumption 
that fih and /Xh are constant on the ring, we obtain 

„2 



7h,i 



mr 



(nl 



/ih) smjcosi. 



(12) 



The potential of an exponential disc in the outer region is 
approximated as [see equation (2P-5) of Binney & Tremain 

(1987)], 



^d{R,z) 



GMd 



3Ri (R^ 



1z' 



2r4 



(13) 



Thus, the y- and z-components of the force due to the inner 
disc per unit mass are, respectively, given by 



and 

Fd,. = ■ 



GA/d y 



GMaz 



1 



1 + 



2z' 



2r^ 



\hR\ [R^ - 2z^) 



2r* 



+ 



37?d 



6Rl 



(14) 



(15) 
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is dominant. Therefore, the precession rate increases with 
decreasing inclination, and vice versa as seen in Fig. m 



fig7.gif 



Figure 7. A sketch of the type I (a) and the type II (b) warps 
seen edge-on. 



Provided that the line of nodes with respect to the equatorial 
plane of the halo is aligned between the inner disc and the 
outer ring, the a;-component of the torque exerted by the 
inner disc is 



Td,x 



mr^ 9RiGMa 



sin 5 cos S, 



(16) 



where S is the inclination angle of the outer ring relative 
to the inner disc plane, and its sign is positive when the 
inclination angle of the inner disc is larger than that of the 
outer ring, and vice versa. 

Next, the absolute value of the total angular momen- 
tum of the outer ring, L, is mr^Q,, where Q is the orbital 
frequency, if the precession rate is rather smaller than the 
orbital frequency. The perpendicular component to the z- 
axis is mr^flsini. The change in Lx over an infinitesimally 
small time At is mr^ Q, sin iujp At, where ujp is the precession 
rate of the outer ring. Thus, 



mr Q sin i ujr, 



(17) 



This should be equal to the torque on the outer ring, so that 
the precession rate is derived as 



Mh 



9RlGMd sin (5 cos 5 



(18) 



'^ 2n 2nr^ sini 

If the warp is of the type I shape such as Fig. fm, i is 
nearly equal to S. It follows from equation (Flq) that the pre- 
cession rate is proportional to cosi. In this case, the second 
term of equation (llq) is negative because 5 is positive, and 
the first term is also negative because Q-t < Hh for the oblate 
halo. Therefore, the precession rate increases with increasing 
inclination as seen in Fig. 0, and the difference in longitude 
between the outer ring and the inner disc becomes larger 
with time. If the warped disc is of the type II shape such 
as Fig. 0b, S is negative. Taking into consideration Qh > ^h 
for the prolate halo, both terms in equation (nm are posi- 
tive. Moreover, i will become smaller than \5\ with decreas- 
ing inclination, so that the second term of equation (|l8[) 



4.3 The inclination 

The time evolution of the inclination is also explained simply 
by the torque on the outer ring. The geometry is the same 
as that used in the previous subsection except that the xy- 
plane is in the precessing frame of the outer ring. 

The 2/-component of the torque, Ty, changes that of the 
angular momentum, Ly — —mr^Qsini, and so, affects the 
inclination angle, i. At the beginning, Ly suffers no change 
from the torque due to the halo because of the axisymmetric 
nature of the halo. If the longitude of the outer ring is the 
same as that of the inner disc, we obtain Ty = 0, and so, the 
inclination remains unchanged. If the longitude of the outer 
ring is smaller than that of the inner disc, i.e., a warp like 
a trailing spiral, Ty becomes positive. Hence, Ly increases, 
so that the inclination angle of the outer ring decreases. 
Since the torque on the inner disc is the opposite direction 
to that of the outer ring according to Newton's third law, 
the inclination of the inner disc increases. As a result, the 
type II warp is generated. Similarly, a warp like a leading 
spiral leads to the type I warp. 

If the precession rate in the prolate halo decreases with 
radius, the warped configuration becomes similar to a trail- 
ing spiral because the sense of the precession is the same di- 
rection as the disc rotation. Consequently, the type II warp 
is produced, which leads to the situation where the incli- 
nation of the outer ring decreases, and the precession rate 
increases. In the subsequent evolution, the longitude of the 
outer ring passes through that of the inner disc. At the point, 
the warped configuration turns into that similar to a lead- 
ing spiral, which, this time, leads to the situation where the 
inclination increases and the precession rate decreases. This 
kind of self-regulation enables the warped disc to avoid the 
differential precession. For the oblate halo, however, once the 
precession of the outer ring recedes from that of the inner 
disc, the difference in precession rate continues to increase, 
so that the warp winds up tightly. 



4.4 Comparison with discrete bending modes 

Our oblate halo model cannot sustain the warped disc, 
though SC found that long-lived warps do exist in some halo 
models. This discrepancy is considered to originate from the 
radial extent of a disc. SC showed that a discrete mode with 
an eigenfrequency u may exist only ii Q — ^ < u) < Q + ^ is 
satisfied at the edge of the disc, where fl and jj, are the or- 
bital and vertical frequencies, respectively. For oblate haloes, 
f7 — /i and u) are negative. Since fi — ^ tends to zero with 
radius, there is some radius where uj is equal to SI — /^. Be- 
yond such a radius, the condition of the existence of dis- 
crete modes is violated. According to SC, bending modes 
become continuous in frequency uj, if a warped disc embed- 
ded in an oblate halo extends beyond the radius at which 
u! = il — fi holds. Consequently, such continuous modes will 
be propagated with a group velocity and disappear (Hunter 
& Toomre 1969). As is found from Fig. ga, such a reso- 
nance radius emerges at _R ~ 4.8 Rd in our oblate halo 
model. This implies that there would be no discrete bend- 
ing mode and that a warped configuration would be coerced 
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1991). These observations appear to favour the view that 
we just witness the different phases of the evolving warped 
discs in prolate haloes, because only warped discs in prolate 
haloes can change the spirality of the line of nodes according 
to our simulations. 

Putting together our simulations and the observations 
mentioned above, we can infer that warps are formed and 
maintained in prolate haloes. However, our fixed halo mod- 
els are quite simplified. In particular, such models cannot 
include the effect of dynamical friction between the warped 
disc and the halo. Dubinski & Kuijken (1995) and Nelson & 
Tremaine (1995) have shown that dynamical friction plays 
an important role to precessing discs. Therefore, we will 
need the simulations of warped discs embedded in live haloes 
to determine the precise evolution of warps, though a huge 
number of particles will be required to incorporate the effect 
of dynamical friction accurately into such simulations, and 
to avoid disc heating due to an insufficient number of halo 
particles. This line of investigation is in progress (Ideta et 
al. in preparation). 



Figure 8. Plots of the frequencies, f2 — /i, f2 + /.i, and oj for 
the oblate halo model (a) and the prolate halo model (b). In the 
calculation of H — ^ and 0-|-/i, the approximated form of equation 
(llSh is used for the disc potential. 



to disperse. Thus, the disappearance of the warping for the 
oblate halo could be due to the existence of the resonance. 
On the other hand, in our prolate halo model, the condi- 
tion ofJl — /x<ci;<r2-|-/iis exactly satisfied within the 
truncated radius, 15 i?d, as shown in Fig. pp. Therefore, a 
discrete mode can exist in our prolate halo model. 

Our adopted disc model is nothing special in the sense 
that the observed light distributions of galactic discs are 
well-described by an exponential law, though a constant 
mass-to-luminosity ratio throughout the disc is assumed. 
Our simulations suggest that as long as a disc is not trun- 
cated abruptly, real galactic discs would not have a discrete 
bending mode if surrounding haloes are oblate. 



5 CONCLUSIONS 

In this paper, we have examined the time evolution of 
warped discs in the oblate and prolate haloes using A'^-body 
simulations. The haloes were represented by fixed external 
potentials in which self-gravitating discs were embedded. 
Then, we have found the warped configurations both in the 
oblate and prolate haloes. While the warping in the oblate 
halo continued to wind up with time and finally disappeared, 
the warping in the prolate halo survived to the end of the 
simulation by regulating the line of nodes of the warped disc 
to be straight. We have shown that this difference in winding 
between the oblate and prolate haloes can be attributed to 
the gravitational torque between the inner and outer discs. 
Observationally, some galaxies show the straight line of 
nodes of the warp within the Holmberg radius, beyond which 
the warp is traced as a leading spiral (Briggs 1990). Others 
show that the line of nodes is delineated as a trailing spi- 
ral (Christodoulou, Tohline, Steiman-Cameron 1988; Bosma 
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